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ABSTRACT 


An  aerodynamic  model  is  constructed  for  the  application 
of  the  properties  of  dynamic  stall  of  airfoils  to  the  calculation 
of  the  airloads  and  blade  motion  of  helicopter  rotor  blades. 
Dynamic  stall  occurs  on  an  airfoil  undergoing  pitc^xing  motion  at 
high  angle  of  attack,  and  is  characterized  by  peak  section  lift 
and  moment  much  larger  than  the  corresponding  static  stall  loads. 
A  method  is  developed  for  the  solution  of  the  equations  of  motion 
of  a  rotor  blade  by  means  of  harmonic  analysis.  The  effect  of 
dynamic  stall  on  the  blade  torsional  motion  at  high  advance  ratio 
is  examined,  and  comparison  is  made  with  the  limited  experimental 
data  available.  An  increase  in  the  dynamic  stall  angle  is  shown 
to  significantly  decrease  the  amplitude  of  the  pitch  motions. 
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NOMENCLATURE 


lift  curve  slope 

blade  sernichord 

section  drag  coefficient 

section  lift  coefficient 

section  moment  coefficient 

Drag 

angle  between  tip  path  plane  and  the  forward  velocity 
positive  for  rearward  tilt  of  the  thrust 

imaginary  number: 

Lift 

Moment 

blade  spanwise  variable 

Rotor  radius 

section  angle  of  attack 

nondimensional  rate  of  change  of  angle  of  attack,  at 
the  instant  of  stall 

blade  flapping  angle 

blade  pitch  angle 

rotor  advance  ratio:  forv;ard  speed  divld«--d  by  QR 

fluid  density 

rotor  azimuthal  variable 

rotor  rotational  speed 
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SECTION  1 


INTRODUCTION 

The  large  transient  control  loads  encountered  at  high 
forward  speeds  or  high  blade  loadings  remain  a  significant  problem 
in  the  operation  and  design  of  helicopters.  A  primary  source  of 
these  loads  is  the  aerodynamic  moment  due  to  blade  stall.  Results 
of  an  investigation  into  the  nature  and  effects  of  stall  in  the 
aerodynamic  environment  characteristic  of  helicopter  blades 
(References  1  to  indicate  that  the  stall  phenomenon  exhibited 
by  an  airfoil  undergoing  a  rapidly  changing  angle  of  attack, 
known  as  dynamic  stall,  differs  substantially  from  the  familiar 
static  stall  at  constant  angle  of  attack.  Such  a  situation  of 
rapidly  Increasing  angle  of  attack  characterizes  a  helicopter 
rotor  blade  a's  it  traverses  the  rotor  disc,  particularly  on  the 
retreating  side  of  the  disc  at  high  advance  ratios.  This  report 
presents  a  method  which  was  developed  to  calculate  the  aerodynamic 
loading  and  resulting  blade  motion  of  a  rotor  blade,  particularly 
the  aerodynamic  moment  and  resulting  blade  torsional  motion  due 
to  dynamic  stall.  The  aerodynamic  model  developed  was  specifically 
designed  for  high  forward  speeds  (y  >  0.25),  where  a  simple  wake 
geometry  model  is  accurate.  Moreover,  due  to  the  lack  of  symmetry 
over  the  rotor  disk  at  high  forward  speed,  recurrent  entry  to  and 
exit  from  stall  (stall  flutter)  is  not  likely  to  occur,  and  the 
very  complicated  and  detailed  wake  and  aerodynamics  due  to  such 
a  motion  need  not  be  considered. 

In  v/hat  follows  all  quantities  are  nondimens ional,  based 
on  the  fluid  density,  the  rotor  radius,  and  the  rotor  rotational 
speed  (p,  R,  fi).  The  section  lift,  moment,  and  drag  coefficients 
are  conventionally  defined. 
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SECTION  2 


THE  AIRLOADS  AND  BLADE  MOTION  CALCULATION 


The  calculation  procedure  used  to  determine  the  airloads 
and  blade  motion  of  a  helicopter  rotor  blade  follows  a  reference 
blade  around  the  azimuth  in  discrete  steps  A'P  of  the  order  of  ten 
degrees.  The  downwash  at  several  blade  stations,  the  section 
loading  at  several  blade  stations,  the  blade  motion,  and  the  wake 
geometry  changes  are  calculated  successively;  then  the  blade  is 
moved  to  the  next  azimuth  station.  This  procedure  is  Illustrated 
in  Figure  1.  The  use  of  a  variable  inflow  calculated  on  the  basis 
of  a  wake  of  trailed  and  shed  vorticity  precludes  any  closed  form 
solution.  The  calculations  are  therefore  performed  iteratively 
until  they  converge  to  the  steady  state  solution.  The  solution 
is  determined  by  the  physical  properties  of  the  blade  (the  semi- 
chord  to  span  ratio  b,  the  Lock  number  L.N.,  and  others),  and 
by  the  root  collective  pitch  0°,  the  tip  path  plane  inclination 
angle  i,  the  advance  ration  y,  and  tLe  tip  Mach  number  ^ 

2.1  The  Aerodynamic  Loading 


The  section  lift,  moment,  and  drag  —  L,  M,  D  —  are 
determined  from  the  appropriate  aerodynamic  theory  based  on  the 
state  of  the  flow  and  the  blade  motion.  The  state  of  the  flow 
(attached  or  nonattached)  at  a  certain  blade  station  is  determined 
by  the  past  history  of  the  flow  and  by  the  angle  of  attack  a: 


^  &  —  Ta/yn 

where 

U-r*  ~  "t  ^ 

lAip  =  \  -h  4  4- 

==  V  -4-  r- 


VA-r 
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(with  the  convention  that  a  positive  angle  of  attack  in  the 
reverse  flow  region  corresponds  to  a  downward  acting  lift).  The 
state  of  the  flow  is  determined  as  follows: 


(a)  For  [aj  less  than  a  defined^  stall  angle,  the  flow  is 
considered  attached.  The  loads  are  determined  using 
lifting  line  and  lifting  surface  theories. 


(b) 


When  |al  becomes  greater  than  some  ct^g,  the  section 
is  considered  to  stall  dynamically.  The  loads  are 
determined  by  the  parameter  ac/V  at  the  instant  of 


stall.  The  dynamic  stall  state  lasts  until  the 
section  angle  of  attack  begins  to  decrease  (but  at 
least  for  an  interval  s  and  then  the  flow  is 
considered  separated.  During  the  separated  state  the 


loads  decay  from  the  dynamic  stall  values  to  static 
stall  values.  When  la|  next  goes  below  some  Ogg  the 
flow  becomes  attached  again. 


(c)  When  |ot(  becomes  greater  than  some  the  flow  is 
called  feathered  (typically  Up  =  60  degrees).  When 
|a(  next  goes  below  Up,  the  flow  is  considered 
statically  stalled,  until  |a|  goes  below  Ugg  and  the 
flow  reattaches. 


The  only  feature  of  the  reverse  flow  region  which  is 
considered  is  the  direction  of  the  component  of  the  velocity  u,p 
normal  to  the  span  direction.  This  determines  the  double  sign: 


which  will  occur  frequently  below. 


Several  additional  points 


must  be  considered: 


(a)  For  points  too  near  the  reverse  flow  boundary  (u^  =  0) 
the  loading  is  set  to  zero. 
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(b)  For  radial  stations  which  at  the  last  azimuth  station 
were  on  the  other  side  of  the  reverse  flow  boundary , 
so  that  the  past  history  of  the  flow  is  not  well 
defined  by  the.,  above  model,  the  flow  is  considered 

to  have  been  attached  at  the  last  azimuth  station. 
Whatever  the  state  of  the  flow  is  near  the  reverse 
flow  boundary,  the  loading  will  be  small  because  of 
the  low  dynamic  pressure  of  the  velocity  normal  to 
the  span. 

(c)  If  the  angle  of  attack  changes  sign  in  the  separated 
region,  so  that  a  more  detailed  calculation  along 
the  azimuth  would  show  the  section  flow  reattaching 
and  then  restalllng  as  the  angle  of  attack  goes 
through  zero,  the  section  is  considered  to  be  stati¬ 
cally  stalled. 

The  logic  described  above  is  shown  in  Figure  2.  This 
model  is  designed  to  show  the  effects  of  a  rather  impulsive  load¬ 
ing  (due  to  dynamic  stall)  on  the  torsional  motion  of  a  blade 
under  high  loading  conditions  and  high  advance  ratio.  It  is  not 
appropriate  for  handling  several  successive  excursions  in  and  out 

of  stall  unless  an  impractlcally  small  azimuth  step  is  used. 

The  convention  for  the  positive  directions  of  the  loads  is 
as  follows: 

lift,  L  upward 

moment,  M  nose  up,  about  the  quarter  chord  (the 

feathering  axis) 

drag,  D  opposing  the  blade  rotation 

pressure,  -Ap  upward 

circulation,  T  directed  outward  along  the  blade 


2.1.1  Attached  Flow 


For  the  blade  section  in  attached  flow,  the  loading  is 
calculated  using  lifting  line  theory;  that  is,  the  flow  is  con¬ 
sidered  locally  two-dimensional,  with  the  influence  of  the  rest 
of  the  blade  and  the  rotor  wake  represented  by  Just  a  downwash 
velocity  at  the  blade  section.  The  aerodynamic  loading  is  due 
to  the  blade  motion  and  the  v;ake  induced  downwash.  The  downwash 
is  calculated  from  the  trailed  and  shed  vorticity  in  the  v;ake  of 
each  blade  of  the  rotor.  In  early  calculations,  only  the  downwash 
due  to  the  trailed  vorticity  was  calculated  in  this  manner;  the 
effect  of  the  shed  vorticity  behind  the  reference  blade  was 
accounted  for  by  the  use  of  Theodorsen's  lift  deficiency  function. 
The  two  methods  of  handling  the  shed  vorticity  showed  little 
difference  as  far  as  stall  loads  and  blade  motion  were  concerned. 

As  a  first  order  correction  for  compressibility  effects, 
the  Prandtl-Glauert  factor  is  used,  based  on  the  Mach  number  of 
the  flow  normal  to  the  span  of  the  blade.  Thus 


L 


incompressible  lifting 
line  theory 


M  =.  -  r  M  \ 

^  1  I  ,  _  J  incompressible  lifting 

'  line  theory 


where 


Mk) 


Mach  number  of  tip  (r  =  r^) 
at  ^  =  90  degrees  ^ 
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2.1.2  Nonattached  Flow 

The  aerodynamics  of  the  flow  when  it  is  stalled,  separatedj 

or  feathered  give  the  section  load  coefficients  —  c, ,  c  ,  c ,  — 

i  m  d 

for  lift  positive  upward,  moment  (about  the  feathering  axis) 
positive  nose  upward,  and  drag  positive  in  the  direction  opposing 
the  relative  velocity.  The  lift  and  drag,  L*  and  D*,  are  respec¬ 
tively  perpendicular  to  and  in  the  direction  of  the  actual  section 
relative  velocity  (see  Figure  3).  Converting  these  coefficients 
to  loads  based  on  and  R  gives 

Ut  ^ 


Referring  again  to  Figure  3,  the  loads  can  be  converted  to  lift 
and  drag  in  the  no  feathering  plane  using 


These  expressions  are  necessary  because  of  the  large  angles  which 
may  be  involved  in  nonattached  flow.  The  circulation  is  given  by 

\  3  ^  LaT  ^ - 

and  the  actual  loads  by  a.^/2-n  times  L  and  M  (for  use  in  the 
equations  for  the  blade  motion  however,  the  lift  curve  slope 
correction  factor  is  already  included  in  the  Lock  number) . 
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2.1.3  Dynamic  Stall 


The  loss  of  lift  and  circulation  and  the  shift  of  center 
of  pressure  an  airfoil  at  an  angle  of  attack  above  -the  static 
stall  angle  (generally  about  8  to  12  degrees)  are  familiar 
aerodynamic  phenomena.  The  operation  of  a  helicopter  rotor  blade 
in  an  environment  of  high  loading  and  changing  angle  of  attack 
makes  it  necessary  to  consider  stall  as  a  dynamic  rather  than  a 
static  phenomenon.  The  unsteady  aerodynamics  involved  in  dynamic 
stall  are  not  yet  completely  understood  (see  References  1  to 
but  fundamental  to  characterizing  the  problem  is  the  time  scale 
of  the  airfoil  motion,  as  represented  by  the  parameter  ac/V.  The 
primary  characteristics  of  dynamic  stall  are  its  occurrence  at  an 
angle  of  attack  greater  than  the  static  stall  angle,  followed  by 
the  shedding  of  vorticity  from  the  leading  and  trailing  edges. 

The  unsteady  aerodynamic  forces  due  to  the  vorticity  passing  over 
the  upper  surface  of  the  airfoil  produce  a  lift  and  nose  dovm 
moment,  with  peak  values  much  greater  than  the  corresponding 
static  stall  loads.  Experimental  and  theoretical  work  (Reference 
has  shown  that  the  primary  parameter  of  dynamic  stall  is  the  rate 
of  change  of  angle  of  attack  at  the  instant  of  stall,  in  the  non- 
dimensional  form 


U- 


with  little  dependence  on  the  actual  angle  of  dynamic  stall. 

Figure  ^  shows  the  results  of  both  theory  and  experiment  for  the 

variation  with  ^^*c/v  of  the  peak  lift  and  moment  coefficients, 

c,  and  c  .  Although  still  requiring  experimental  verlfica- 
■*’max  max 

tion,  the  theory  of  Reference  implies  that  the  peak  coefficients 

approach  limits  of  the  order  of  3.0  for  c,  and  -0.8  for  c 

i  m 

max  max 


for  high  stall  rates. 
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As  a  model  for  dynamic  stal3 ,  It  is  assumed  that  above  a 
given  angle  of  attack,  as  long  as  the  angle  of  attack  is 

increasing,  the  lift  and  moment  coefficients  are  maintained  at 
constant  values.  These  coefficients  are  taken  as  equal  to  the 
peak  coefficients  as  given  by  the  experimental  and  theoretical 
results  cited  aobve,  and  are  therefore  dependent  on  the  stall 
rate  parameter  oeH/V.  The  empirical  representation  used  is  also 
shown  in  Figure  4.  Data  (Reference  4)  indicate  that  the  dynamic 
stall  angle  is  about  14  to  20  degrees;  typically  the  value  of 
“ds  “  degrees  is  used.  Thus  the  value  of  at  the  instant 

of  stall  gives 

Cw,,= 

The  use  of  this  model  (Reference  5)  has  shown  that  typically 
the  large  nose  down  moment  generated  by  dynamic  stall  quickly 
causes  the  angle  of  attack  to  begin  to  decrease  (due  to  elastic 
pitch  motion  of  the  blade)  and  thus  the  flow  to  enter  the  separa¬ 
tion  model.  The  dynamic  stall  region  extends  only  for  about 
15  degrees.  The  resultant  elastic  pitch  excitation  is 
characteristically  very  impulsive.  At  high  advance  ratios  the 
model  is  only  weakly  dependent  on  the  stall  angle  a^g.  The 
possibility  of  self-excited  elastic  torsional  motions  (stall 
flutter)  is  likely  only  at  very  low  advance  ratio,  with  a  great 
deal  of  symmetry  around  the  rotor  disk.  At  very  high  advance 
ratios  (y  >  0.5),  the  value  of  is  frequently  larger  than 

the  range  shown  in  Figure  4,  and  so  the  dynamic  stall  loads  are 
very  large. 

The  use  of  harmonic  analysis  to  represent  the  load  data, 
which  is  calculated  only  at  discrete  azimuth  points,  as  a  smooth 
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distribution  around  the  disk  actually  replaces  the  instantaneous 
rise  of  the  stall  loads  to  the  peak  dynamic  stall  values  by  a 
gradual  increase,  with  the  rise  time  determined  by  the  azimuth 
increment  in  the  harmonic  analysis  (see  Section  2.5.2).  For  a 
very  small  azimuth  increment  in  the  calculation  process,  some  limit 
must  be  placed  on  the  rate  of  increase  of  the  stall  loads.  Using 
a  linear  rise.,  the  lift  and  moment  coefficients  are 


for  <r'  ^ 


-h 


ys.KVg  , 

<^vv\  ~  ^ -f*  c 

and  for 

—I  C■^A^  Q 


where  c,  and  c  are  the  coefficients  for  the  attached  flow  at 
a  a 

the  dynamic  stall  angle  and  is  the  azimuth  increment  since  the 
dynamic  stall  angle  was  reached.  The  experimental  and  theoretical 
data  (Reference  4)  suggest  that  the  rise  time  should  be 
approximately  10  degrees,  which  is  also  confirmed  by  a  comparison 
of  the  results  of  this  calculation  model  and  experimental  data. 


2.1.^  Separation 


Experimental  results  (Reference  6)  show  that  the  flow  over 
an  airfoil  oscillating  at  high  angle  of  attack  separates  at  the 
peak  angle  of  attack  for  values  of  the  reduced  frequency,  k  =  wb/V, 
greater  than  0.3.  It  is  also  shown  there  that  the  Theodorsen 
theory  of  oscillating  airfoils  is  equivalent  to  a  first  order 
differential  equation  for  the  moment  coefficient  of  the  form 


V'A— 


C  C  VM  - C'vw.  — 


\ 

-T^ 
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where  the  subscript  i  refers  to  the  Initial  conditions,  and  the 
Theodorsen  theory  gives  X  and  the  time  constant  as  functions 
of  k  and  C(k)  (the  lift  deficiency  function).  A  similar  expression 
may  be  derived  for  the  lift  coefficient,  with  a  time  constant  T^. 
When  this  equation  was  assumed  to  hold  in  the  stalled  region  also, 
moment  variation  with  angle  of  attack  was  calculated  which  was 
in  go  3d  agreement  with  the  experimental  results.  Since  a  first 
order  equation  would  be  expected  to  give  a  reasonable  approxima¬ 
tion  to  the  decay  of  the  dynamic  stall  loads  (see  Reference  3  for 
experimental  data  for  this  decay)  this  agreement  primarily  indicates 
that  the  time  constants  given  by  the  Theodorsen  theory  may  also 
be  used  when  the  airfoil  is  stalled. 


Therefore  as  a  model  for  separated  flow,  it  is  assumed  that 
after  the  section  angle  of  attack  reaches  its  maximum  value,  the 
lift  and  moment  coefficients  decay  exponentially  to  static  stall 
values.  Thus  the  lift  and  moment  coefficients  are 


+  CS5 


C.vv\  —  <2^ 


4'  C 


Where  c^  and  c  are  the  coefficients  at  the  peak  angle  of  attack 
o  o 

(dynamic  stall  coefficients)  and  c,  and  c  are  the  static  stall 

s  s 

coefficients  (based  on  the  same  sign(a)  and  sign(Um)  as  are  c, 

■^o 

and  c  );  %  is  the  azimuth  angle  at  which  the  separation  region 
^0 

was  entered.  The  time  constants  are  given  by  the  Theodorsen  theory 
as  approximately  =  1.0  and  =  2.5.  Calculations  with  these 
values  show  that  the  time  constants  used  for  the  transient  stall 
of  the  rotor  blade  should  be  several  times  these  values,  which 
were  derived  and  experimentally  confirmed  for  a  highly  periodic 
stalling  motion. 
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Dynamic  stall  as  described  above  is  a  phenomenon  rather 
.'ndependent  of  the  airfoil  motion  after  its  initiation  at  the 
dynamic  stall  angle.  Thus  the  dynamic  stall  state  should  exist 
for  a  definite  period,  regardless  of  the  subsequent  blade  motion. 
Therefore  the  separation  state  is  assumed  not  to  start  until  at 
least  after  the  dynamic  stall  angle  is  reached;  if  the  angle 
of  attack  is  still  increasing  after  this  point,  the  separation 
state  begins  when  the  angle  of  attack  starts  to  decrease. 

It  is  assumed  that  the  flow  reattaches  at  the  static  stall 
angle  a^g.  However,  reattachment  is  also  delayed  if  necessary 
until  after  the  initiation  of  dynamic  stall. 

2.1.5  Static  Stall 

The  model  used  for  static  stall  assumes  that  the  lift  and 

moment  coefficients  are  maintained  at  the  constant  static  stall 

values,  c,  and  c  .  Then  the  coefficients  are 
s  s 

ojj  =  .=!=.  (c<) 

OvvA  =  •4- 

Typical  values  used  are  c,  =1.0  and  c„  =  -0.15,  and  for  the 

S  5 

static  stall  angle  Ugg  =  12  degrees. 

2.1.6  Feathered  Flow 

The  experimental  results  of  Reference  7  shov/  that  the  stall 
model  is  reasonably  accurate  up  to  angles  of  attack  of  around  6o 
degrees.  These  experiments  were  static  measurements  on  a  NACA 
0012  airfoil  section,  A  good  approximation  to  the  experimental 
coefficients  above  the  feathered  angle  is  given  by 

—  i  oCi 
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for  la|  >  a_.  The  experiments  give  the  parameters:  c,  =  c,  =1.0; 

^  "^S 

=  -0.40;  and  Op  =  60  degrees. 

2.1.7  The  Drag 

The  experimental  data  of  Reference  7  also  give  the  drag 
coefficient  for  the  NACA  0012  airfoil  section.  Unlike  for  attached 
flow,  for  nonattached  flow  the  drag  may  have  a  significant  component 
in  the  rotor  thrust  direction  because  of  the  large  local  angle  of 
attack.  For  nonattached  flow,  the  data  is  approximated  quite  well 
by:  ^ 

<=-A  —  2,*0  (  S'.M  c<'> 


2.2  The  Rotor  Wake  and  Downwash 

The  wake  of  a  helicopter  rotor  blade  in  forward  flight 
consists  of  trailed  and  shed  vortlcity  in  a  distorted,  skewed, 
helix.  In  order  to  calculate  the  dovmwash  induced  at  the  rotor 
blade  by  this  v/ake,  it  is  represented  by  a  net  of  finite  strength, 
finite  length  line  vortices  behind  each  blade.  The  distribution 
of  the  blade  bound  circulation  over  the  rotor  disk,  r(r,''lf^),  is 
sufficient  to  establish,  by  continuity  of  vortex  lines,  the 
strength  of  each  line  element  in  the  net.  A  semi-rigid  wake 
geometry  model  is  used  to  give  the  position  of  each  line  element; 
that  is,  each  point  of  the  wake  is  assumed  to  travel  downward 
always  at  a  velocity  equal  to  the  downwash  at  the  point  on  the 
disk  from  which  it  was  trsiiled  or  shed.  Prom  its  strength  and 
position,  the  downv;ash  Induced  by  each  line  element  may  be  calcu¬ 
lated;  the  downwash  at  the  rotor  blade  is  then  the  sum  of  the 
contributions  from  all  elements  in  the  wake. 

2.3  The  Reverse  Flow  Region 

It  is  shown  in  Reference  7  that  an  airfoil  in  reverse  flow 
has  almost  the  same  lift  and  moment  variation  with  angle  of  attack 
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as  for  conventional  flow,  as  would  be  expected  with  .hin  airfoils 
Therefore  the  only  change  in  the  model  is  in  the  positive  direc¬ 
tions  of  the  loads.  As  far  as  the  blade  torsional  motion  is 
concerned,  the  most  important  feature  of  the  reverse  flow  region 
is  that  the  aerodynamic  center  is  then  at  the  three-quarter  chord 
so  the  lift  gives  a  large  moment  about  the  blade  feathering  axis, 
normally  located  at  the  quarter  chord. 

2.^  Compressibility  effects 

The  only  consideration  of  compressibility  effects  on  the 
loads  is  the  use  of  the  Prandtl-Glauert  correction  factor  based 
on  the  Mach  number  of  the  local  velocity  perpendicular  to  the 
span,  u,p.  The  neglect  of  the  effect  of  compressibility  on  the 
stall  loads  is  a  serious  defect,  as  it  is  well  established  that 
this  is  a  significant  factor  for  airfoils  stalling  in  the  environ 
ment  characteristic  of  helicopter  blades.  However,  once  experi¬ 
mental  and  theoretical  data  on  the  effects  of  compressibility  on 
dynamic  stall  become  available,  it  should  be  easy  to  incorporate 
them  into  the  aerodynamic  model  which  has  been  described  here. 

2.5  The  Blade  Motion 

The  vertical  position  of  the  blade  is  composed  of  the 
feathering  axis  deflection,  z^,  and  a  chordwise  rigid  pitch 
motion,  0: 

~  -I-  X 

The  deflection  of  the  elastic  axis  is  represented  as  th  sum  of 
rigid  flapping  and  first  mode  elastic  bending;  and  'che  pitch 
motion  as  tre  sum  of  a  rigid  (spanwise)  pitching  mode  and  first 
mode  elastic  pitching: 

e  ecf 
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where  ?(0)  =  0.  An  approximate  mode  shape  is  used  for  elastic 
bending,  and  the  elastic  mode  is  typically  linear  twist: 


v^Cr')  ^  p.  2-  „  3  ^ 
;=r  V- 


The  equation  of  motion  is  most  conveniently  solved  in  the  no 
feathering  plane,  where  the  only  control  input  is  collective 
pitch  (the  tip  path  inclination  angle  i  also  enters  indirectly 
through  the  downwash).  Cyclic  pitch  can  be  varied  arbitrarily 
with  the  orientation  of  the  shaft;  the  dynamics  and  aerodynamics 
of  the  rotor  determine  the  orientation  of  the  no  feathering  plane 
with  respect  to  the  tip  path  plane. 


2,5,1  The  Equations  of  Motion 

The  equations  of  motion  for  a  constant  chord  (uniform 
spanwise  structural  properties),  coincident  elastic  axis  and 
inertial  axis,  articulated  (no  lag  hinge,  zero  flapping  hinge 
offset)  blade  may  be  derived  following  Reference  8: 

=  L.N. 


(T  ©  “f-  ©  4-  Co©  ^  ©  —  ©*  4'  C  © 


r«r 

-  L.h.  V  . 


M 


4-  © 


T4 


1*1 


where 


The  lift  and  moment,  L  and  M,  are  based  on  the  ideal  lift  curve 
slope  2tt;  the  Lock  number  includes  the  lift  curve  slope  correc¬ 
tion  factor. 

The  pitching  equations  of  motion  are  very  sensitive  to 
impulsive  moment  excitation,  and  in  order  to  obtain  a  convergent 
solution  by  means  of  harmonic  analysis,  it  is  necessary  to  include 
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part  of  the  aerodynamic  damping  on  the  left  hand  side  of  the 
equations.  Thus  the  nonharmonic  (no  dependence  on^),  incom¬ 
pressible  damping  is  extracted  from  the  lifting  line  result  for 
the  aerodynamic  moment.  To  the  total  (compressible)  moment, 
M/^TTb,  is  added  the  quantity: 


nr 


-4TTr)o 


\  ^ 

r"  ftr  b 


© 


for  the  equations  of  motion  for  B  and  respectively.  Corres¬ 

pondingly  to  the  left  hand  side  of  the  equations  are  added  the 
terms: 


icpe 


and 


where 


(The  damping  terms  given  here  are  for  the  feathering  axis  at  the 
aerodynamic  center  only.) 

2.5.2  Solution  for  the  Blade  Motion 

The  equations  of  motion  are  most  conveniently  solved  by 
harmonic  analysis  of  the  loads  and  blade  motion.  Thus  writing 
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and  similarly  decompcslng  Qj  6,  0^,  8°,  0^,  and  the  integrals  of 
the  lift  and  moment,  the  equations  of  motion  become: 


=  L.k>. 

(:  <?r  -  =  L.K>. 


0  2.,/  M 

+  L.f^. 


=  e“  nz 

vk 


XI:,  ^ 


All  but  the  stead.v  components  (n  =  0)  of  the  collective  and  twist 
inputs,  0°  and  0^,  are  usually  zero.  These  equations  are  easily 
solved  for  the  harmonics  of  the  blade  motion. 


The  right  hand  sides  of  these  equations  are  the  harmonics 
of  the  integrals  over  the  span  of  the  aerodynamic  forcing  loads. 
In  the  calculation  of  rotary  wing  airloads  these  Integrals  will 
be  evaluated  only  at  a  finite  number  of  azimuth  stations  around 
the  disk.  Calling  any  of  these  Integrals  F(^) ,  its  value  is 
calculated  at  J  azimuth  stations.  Then  the  harmonics  of  F  are: 


i- 

■S' 


£ 

a>y.vwUVV, 


r.  d 


t) 


AV.  =  i  ^ 

'a  ^  s 
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so  that 


For  N  =  J/2  this  gives  exactly  P(^j)  =  Pj-j  for  N  <  J/2  this  is  an 
approximation  in  the  least  squares  sense. 

Because  of  the  smoothing  nature  of  the  harmonic  representa¬ 
tion,  the  rate  of  change  of  a  function  between  two  calculated 
points  will  be  several  times  the  rate  of  change  of  a  straight 
line  representation  between  the  two  points.  This  can  introduce 
convergence  problems  when  large  changes  in  the  loads  are  involved, 
as  when  stall  occurs,  because  the  effective  rise  time  of  the  loads 
js  increased.  Therefore  all  the  harmonics,  F^,  are  multiplied  by 
the  factor: 

r  ^  -|2 

J 

The  harmonic  representation  is  thereby  converted  to  a  function 
connecting  the  calculated  points,  Pj  at  by  straight  lines. 

The  exact  representation  requires  an  infinite  series,  so  truncat¬ 
ing  the  series  has  a  smoothing  effect.  More  importantly,  it  can 
be  seen  that  this  factor  decreases  the  magnitude  of  the  harmonics 
of  the  least  squares  representation.  Effectively,  the  rise  time 
of  the  function  is  now  restricted  to  at  least  the  azimuth  incre¬ 
ment,  2it/J,  in  the  harmonic  representation. 

The  use  of  harmonic  analysis  gives  a  smooth  distribution 
of  the  airloads,  which  are  actually  calculated  only  at.  discrete 
points  around  the  azimuth.  This  smoothed  loading  distribution 
then  gives  the  blade  motion  through  the  equations  of  motion.  The 
use  of  harmonic  analysis  has  an  important  consequence  in  the  stall 
model  however.  In  the  stall  model  described  above,  the  lift  and 
moment  coefficients  Jump  immediately  to  the  stall  values  when  che 
section  angle  of  attack  exceeds  the  dynamic  stall  angle.  The  use 
of  harmonic  analysis  smooths  this  Jump  over  the  azimuth  increment. 
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2tt/J,  which  becomes  a  lower  bound  on  the  rise  time  of  the  stall 
loads.  With  a  small  azimuth  increment,  it  is  necessary  to  include 
in  the  stall  model  itself  a  restriction  on  the  rise  time  of  the 
stall  loads  (as  in  Section  2.1.3). 


2.5.3  The  Solution  for  First  Harmonic  Flapplnp; 

Since  the  natural  frequency  of  an  articulated  blade  with 
zero  flapping  hinge  offset  is  exactly  1/REV,  the  equation  of 
rigid  flapping  motion  is  singular  for  the  solution  for  3-j^.  For 
n  =  1  the  equation,  with  the  integral  replaced  by  a  sum  over  the 
blade  span,  reduces  to: 


which  is  to  be  solved  for  3^.  The  values  of  L  over  the  disk  are 

known,  calculated  at  each  azimuth  station  using  some  32j  and  3^^ . 

The  solution  will  be  found  for  3^^  assuming  that  the  contributions 

to  L  from  all  factors  except  those  involving  3^j  and  3qj  will  be 

the  same  for  the  new  values  of  3^^  and  3^.  That  is,  factors  such 

as  changes  in  the  loads  due  to  the  changes  in  the  downwash  associated 

with  the  new  values  of  3,  and  3^  are  considered  as  second  order 

i  o 

effects.  These  and  similar  effects,  such  as  changes  in  stall 
region  boundaries,  and  changes  in  the  pitch  motion,  will  show  up 
eventually  of  course,  so  that  the  solution  for  a  new  3^  and  3^ 
must  be  obtained  at  each  azimuth  station  as  the  calculation 
procedure  goes  around  the  disk,  until  they  converge  to  the  correct 
values. 


These  contributions  to  L  that  are  assumed  to  be  independent 
of  32  and  3^  are  contained  in  the  quantity 


r 


V- 
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where  the  subscript  ’’present"  means  based  on  the  old  values  of  6, 
that  is  on  and  ,  This  is  the  quantity  usually  used  on  the 
right  hand  sides  of  the  equations  of  motion  in  the  solution  for 
the  harmonies  of  the  blade  motion.  Prom  this  must  be  subtracted 
the  contributions  due  to  B-j^j  and  B^j .  The  direct  dependence  of 
L  on  B  enters  only  through  the  lifting  line  theory  results,  and 
since  this  dependence  is  linear,  may  write 

where  St*  !tf.  ,  j  y  are  functions  of  r  and  ^9  »  ip,  ,  ip, 

are  given  by  the  lifting  line  theory  and  includes  all  other 
loads. 


where  as  before  the  subscript  "j"  refers  to  the  sum  over  the 
azimuth  (which  enters  in  obtaining  the  first  harmonic  of 

(  r(L/4‘rrb)Ar) ) .  The  solution  of  this  equation  is: 

spc^tn 

Sl+Xw>S^^+ 
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where 


An  -.\A/\  — 


.«:^^VWXj<-\^ 


(fe^ 


~^^oCSL-iSL) 

^  ^  '  ^*'**'^ 


A.  2 
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and 

-Si— •3=^„-h.  ^  M'j 

sa-  c^s 

,£L 

5^  ^“.^'-p''  s^  '^s  M'lj  f5u>-S(iJp|^ 
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S^  '3'  ..^«t(.  cos  X 
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is  given  as  usual  by 


and  the  solution  for  3, 


U 


Av^ 

>^0  »v 


Now  5^^  and  5?^,  are  given  by  lifting  line  theory  if  the 
flow  is  attached.  Extracting  the  terms  proportional  to  3^,  3-;L  > 
and  3q,  including,  the  component  of  the  downwash  due  to 

the  tilt  of  the  tip  path  plane  with  respect  to  the  no  feathering 
plane,  gives: 

::F  ^  Ut  ('/Acos^— lo')  "] 

*-f“  "’Tf"  W  -f.  ^ 
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•  t)  ^  ^  r 

^  ^  ^  L/ 
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if  the  flow  is  attached,  and 

-5^0  --Sp,  —  o 

if  the  flow  is  nonattached. 

With  the  straight  line  rather  than  the  least  squares 
representation  in  the  harmonic  analysis,  it  is  necessary  to  replace 
the  factor  l/j  by  the  factor: 


where  it  appears  in  the  above. 
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SECTION  3 


RESULTS  . 

3.1  Effect  of  Dynamic  Stall  on  the  Blade  Torsional  Motion 

Reference  5  reports  the  results  of  using  the  aerodynamic 
model  described  in  Section  2  to  calculate  the  airloading  and 
blade  torsional  of  a  helicopter  rotor  at  high  advance  ratio. 

These  results  may  be  summarized  as  follows. 

For  0.^  <  u  <  0.8,  a  stall  region  appears  on  the  retreating 
side  of  the  disk,  outside  the  reverse  flow  region.  The  blade 
pitch  motion  is  characterized  by  an  initial  nose-dovm  impulse  due 
to  dynamic  stall,  a  decrease  in  pitch  in  tne  fourth  quadrant  as 
increasing  dynamic  pressure  increases  the  moment  due  to  separated 
flow,  and  dynamic  overshoot  and  damped  oscillations  in  the  first 
quadrant  after  flow  reattachment  occurs.  The  pitch  oscillations 
are  well  damped  by  potential  flow  aerodynamic  moments.  For  forward 
disc  inclination,  the  stall  area  appears  initially  at  outboard 
blade  stations,  resulting  in  larger  stall  moments  and  pitch  motions. 

For  y  >  0.8,  the  blade  stalls  before  entering  and  after 
leaving  the  reverse  flow  region,  which  occupies  a  large  part  of 
the  retreating  side  of  the  disk.  For  an  articulated  blade,  stall 
is  not  an  Important  feature  of  the  loads  at  extremely  high  advance 
ratio,  because  the  large  collective  pitch  necessary  to  get  signifi¬ 
cant  stall  regions  results  in  unusably  small  or  even  negative 
rotor  thrust. 

For  y  >  0.6,  the  blade  pitch  motion  is  characterized  by  a 
large  pitch  increase  and  decrease  centered  about  ^  =  270  degrees, 
caused  by  the  downward  lift  force  in  the  reverse  flow  region 
acting  at  the  three-quarter  chord. 

The  Importance  of  the  concept  of  dynamic  stall  for  a  rotor 
operating  at  high  advance  ratio  lies  in  the  fact  that  the  rotor 
blade  stalls  at  high  values  of  the  parameter  »  resulting 
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in  lift  and  moment  initially  several  times  larger  than  the 
corresponding  static  stall  values.  It  should  be  noted  also  that 
the  dynamic  overshoot  due  to  flow  reattachment  following  stall 
often  is  the  greatest  blade  motion  due  to  stall,  and  the  severity 
of  these  oscillations  is  likely  to  be  lessened  only  by  an  increase 
in  pitch  damping,  possibly  by  mechanical  means. 

A  presentation  of  the  airloads  and  pitch  motion  for  a  rotor 
in  various  operating  conditions  may  be  found  in  Reference  5. 

3.2  A  Comparison  with  Experiment 

Reference  9  presents  limited  experimental  data  on  the 
torsional  motion  of  a  rotor  blade  in  several  operating  conditions. 
The  theoretical  torsional  motion  was  calculated  for  these  cases 
using  the  model  described  in  Section  2. 

Figure  5  shows  the  pitch  motion  for  y  =  0.68,  i  =  0,  and 
0°  =  iJ.O  degrees  with  the  blade  center  of  gravity  at  the  30  percent 
chord.  This  is  a  ease  of  classical  flutter  due  to  center  of 
gravity-elastic  axis  offset.  Figure  6  shows  the  pitch  motion  for 
]i  -  1.^7,  i  =  0,  and  0°  =  2.0  degrees.  This  is  a  case  with  large 
pitch  motion  due  to  the  lift  in  the  reverse  flow  region.  For 
both  of  these  cases,  no  elastic  pitch  motion  was  allowed. 

Figure  7  shows  the  pitch  motion  for  y  =  C.29^,  i  =  0,  and 
0°  =  11.0  degrees.  This  is  a  case  involving  dynamic  stall  of  the 
blade  on  the  retreating  side  of  the  disk.  Two  theoretical  results 
are  shown:  one  for  a^g  =  1^  degrees  and  no  elastic  pitch  motion; 
and  one  for  =  16  degrees,  delayed  by  20  degrees  (the 

abscissa  for  this  line  is  actually  ^+20°),  with  no  rigid  pitch 
motion  and  an  elastic  pitch  mode  shape  ?  =  sln^r**  . 

A  discussion  of  these  results  may  be  found  in  Reference  10. 
V/ith  so  little  experimental  data  it  is  not  actually  possible  to 
make  an  adequate  evaluation  of  the  dynamic  stall  model  used  here. 
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3.3  The  Effect  of  the  Dynamic  Stall  Angle 

The  dynamic  stall  concept  suggests  a  means  to  lessen  the 
effects  of  stall.  An  increase  in  the  dynamic  stall  angle  will 
result  in  a  significant  reduction  in  the  severity  of  the  pitch 
motions  due  to  dynamic  stall;  indeed,  for  a  given  collective 
pitch,  a  large  enough  dynamic  stall  angle  could  prevent  the 
occurrence  of  stall  at  all.  The  dynamic  stall  angle  is  dependent 
on  the  parameter  airfoil  heaving  motion,  and  on 

airfoil  shape  (References  3  and  ^),  and  thus  airfoils  may  be 
designed  with  the  aim  of  minimizing  stall  effects. 

In  order  to  evaluate  the  effect  of  increasing  the  dynamic 
stall  angle,  the  blade  pitch  motion  was  calculated  for  the  same 
rotor  used  in  the  calculations  of  Reference  5,  with  several  values 
of  Cpg.  The  natural  pitch  frequency,  w^,  was  5.5/REV.  Figures  8, 
9,  and  10  show  the  harmonics  of  the  rigid  pitch  motion,  0^  (for 
simplicity,  no  elastic  pitch  motion  was  considered)  for  three 
operating  conditions:  Figure  8  is  for  y  =  0.5,  i  =  0,  and  0°  =  12 
degrees;  Figure  9  is  for  y  =  0.5,  i  =  0,  and  0°  =  15  degrees;  and 
Figure  10  is  for  y  =  0.3>  1=0,  and  0®  =  11  degrees.  Figure  11 
shows  the  actual  pitch  motion  for  the  same  case  as  in  Figure  10. 
The  decrease  in  the  amplitude  of  the  pitch  motion  as  the  dynamic 
stall  angle  is  increased  is  quite  substantial. 
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SECTION  i} 


CONCLUSION 

An  aerodynamic  and  dynamic  model  has  been  constructed  which 
has  proved  useful  in  calculating  the  torsional  motion  of  a  heli¬ 
copter  rotor  blade  under  operating  conditions  involving  blade 
stall.  The  concept  of  dynamic  stall  and  the  use  of  the  experi¬ 
mental  and  theoretical  data  on  dynamic  stall  are  important  parts 
of  the  aerodynamic  model  used  here.  It  is  clear  that  further  work 
should  be  done  to  determine  the  nature  of  dynamic  stall,  in  order 
that  a  refined,  and  therefore  more  widely  applicable,  model  of 
stall  may  be  constructed;  it  should  be  relatively  simple  to 
incorporate  new  data  into  the  calculation  procedure  described 
here.  There  is  also  a  great  need  for  more  experimental  data  on 
the  torsional  motion  of  rotor  blades,  in  order  that  the  model 
constructed  may  be  rationally  evaluated.  Calculations  of  the 
harmonics  of  the  pitch  motion  for  varying  dynamic  stall  Angle 
have  shown  that  the  amplitude  of  the  torsional  motion  is  significantly 
decreased  by  an  increase  of  the  stall  angle.  This  illustrates  the 
v^lue  of  designing  airfoils  for  good  stall  characteristics. 
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FIGURE  5  BLADE  PITCH  MOTION:  ^=0.5,  i  =  O,  4,  AND  BLADE  CENTER  OF 
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FIGURE  6  BLADE  PITCH  MOTION 
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An  aerodynamic  model  is  constructed  for  the  application  of 
the  properties  of  dynamic  stall  of  airfoils  to  the  calculation  of  the 
airloads  and  blade  motion  of  helicopter  rotor-  blades.  Dynamic  stall 
occurs  on  an  airfoil  unaergoing  pitching  motion  at  high  angle  of 
attack,  and  is  characterized  by  -•eak  section  lift  and  moment  much 
larger  than  the  corresponding  static  stall  loads.  A  method  is 
developc  '  for  the  solution  of  the  equations  of  motion  of  a  rotor 
blade  by  means  of  harmonic  analysis.  The  effect  of  dynamic  stall  on 
the  blade  torsional  motion  at  high  advance  ratio  is  examined,  and 
comparison  is  made  with  tb*'  limited  experimental  data  available.  An 
Increase  in  the  dynamic  st. ''1  angle  is  shown  to  significantly  decrease 
the  amplitude  of  the  pitch  motions. 
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